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1 ABSTRACT / SUMMARY 


This document 

presents 

the source 

code for 

the latest production 

version. 

MARK 

IV, of 

the NASA 

- Langley 

Vortex Lattice Computer 

Program. 

All 

viable 

subcritical 

aerodynamic 

features of previous 

versions 

have 

been retained. This version 

extends the previously 


documented program 

capabilites 

to four planforms, 400 

panels, and 

enables the user 

to obtain 

vortex-f low aerodynamics 

on 

cambered 

planforms, flow field properties 

off the configuration 

in 

attached 

flow, and planform 

longitudinal 

load distributions. 




2 INTRODUCTION 


The NASA - Langley Vortex Lattice FORTRAN Program (VLM) is designed to 
estimate the subsonic aerodynamic characteristics of up to four complex 
planforms. The concepts embodied in this program are mostly detailed 
in references 1,2 and 3; this document is intended to serve as an 
update to these references for users and computer specialists who have 
an interest in implementing this program on their computers. 
Basically, the VLM Program is a segmented program designed to run on 
the Control Data Corporation (CDC) computers with the NOS operating 
system. This program requires a run-time field length of approximately 
130K (octal) words of memory, and uses the Langley Research Center 
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Graphics Output System in FORTRAN (LRCGOSF) , along with numerous CDC 
routines that provide random access of mass storage files. The User s 
Guide associated with this program is listed as reference 4. 


Use of trade names or names of manufacturers in this report does not 
constitute an official endorsement of such products or manufacturers, 
either expressed or implied, by the National Aeronautics and Space 


Administration. 
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3 PROGRAM STRUCTURE 


3.1 OVERVIEW 


The VLM program consists of one FORTRAN program, 39 FORTRAN 
subroutines, seven labeled common blocks and nine files. The overall 
structure of the code can be represented as a tree diagram with seven 
major branches, as follows: 


(2) 

(3) 

(4) 

(5) 

(6) 

AERODYNAMIC 

AERODYNAMIC 

NEAR-FIELD 

SIDE-EDGE 

VORTEX 

EQUATIONS 

ANALYSIS 

SOLUTION 

SUCTION 

FLOW 

MATRIX 

• 

• 

ANALYSIS 

ANALYSIS 

SOLVER . 

• 

• 

• 

• 

# 

• 

# 

• 

• 

(i) 

• 

• • 

• 

• 

• • 
• • 

(7) 

GEOMETRIC 

• • 

• 

• • 

LONGITUDINAL 

COMPUTATIONS 

• • 

» 

• • 

LOAD 

• 

• • 

• ♦ 

• 

ANALYSIS 

• 

« 

• • • 

• 

• 

• if 

• • » • 

• • » • • 

• • • 

• 

• 

• 

PROGRAM 

• 

• • • • • 

• 

• 

» 


WINGAL 
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The main program, WINGAL, forms the "root node" to the tree. Each of 
the seven major branches, or nodes, has a specific function in the 
overall computational process, and some of these branches consist of 
more than one routine (this will be explained in the next section). 
The point of this tree diagram is that it constitutes the way the 
program is loaded into memory by the operating system, and hence, the 
way it executes for any given input set. The actual file of loader 
directives file based on this structure will be discussed in Section 4. 

3.2 FUNCTION BY NODE 


The root node consists of the main program (WINGAL) and four 
subroutines. The main program is used to declare all files and common 
blocks, as well as direct the overall processing done by this code. 
The four subroutines contained in this node are LOADING, FTLUP, INFSUB 
and READIN. LOADING is not called by any routine in VLM; its sole 
purpose is to force the loading of the CDC mass storage routines at 
this level so they are accessible to all higher level nodes that 
require them. FTLUP and INFSUB, however, are used by several of the 
higher nodes, and hence, must be in memory at all times. FTLUP is a 
linear interpolation routine that uses the CDC Fortran routine LOCF to 

determine the absolute location (address) of a variable. READIN is 
used to read in and print out the input data, with line numbers, for 
the user's reference. The common blocks declared in WINGAL serve as 
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the principal method of information transfer between the higher nodes; 
they are considered "global” and must also remain in memory at all 
times. The higher nodes in the tree are moved in and out of central 
memory by the loader when they are called by the main program, and this 
overlapping results in considerable reduction of central memory 
requirements. The function and associated routines of the higher nodes 
are as follows: 


3.2.1 GEOMETRIC COMPUTATIONS 


This node consists of two subroutines, GEOMETRY and PLNPLT. Subroutine 
GEOMETRY is called by WINGAL to determine, from the input data, the 
geometry of the configuration. GEOMETRY then makes a call to PLNPLT, 
which, in turn, produces the "printer plot" of the configuration. 


3.2.2 AERODYNAMIC EQUATIONS MATRIX SOLVER 


This node consists of six routines and performs the Given's Method of 
Matrix Solution to determine the circulation terras of the horseshoe 
vortices; i.e., these routines solve the matrix of the basic linear 
aerodynamic equations. The routines in this node are: 

1. MATXSOL- Called by WINGAL, MATXSOL is the main routine in this 
node and generates the elements in the aerodynamic influence 
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coefficient matrix. MATXSOL then calls routine GIVENS to 
effect a solution. 

2. GIVENS - Called by MATXSOL to partition the work storage 
arrays into rows and columns. 

3. BLOCKR - Called by GIVENS to compute the size of, and the 
number of rows in each partition of the triangularized matrix. 

4. TRIANG - Called by GIVENS to triangularize the augmented 
matrix using planar rotations. 

5* SOLVER - Called by TRIANG to perform the back substitution on 
the matrix and overstore the results back onto mass storage. 

6. BUFFIN - Called by TRIANG to transfer data into the work array 
from mass storage. 


3.2.3 AERODYNAMIC ANALYSIS 


This node consists of five routines and is used to compute the linear 
aerodynamic characteristics of the configuration. These routines are: 


1. AERODYN - Called by WINGAL, AERODYN is the main routine in 
this node. AERODYN computes the linear aerodynamic lifting 
pressures and overall forces and moments. 

2. FLOWFL - Determines the flow field characteristics off the 
wing. FLOWFL reads the field line definition when the flow 
field data is required. FLOWFL also uses routine FTLUP. 

3. HEAPSRT - Called by FLOWFL to sort the X-Y values of the panel 
centers. FLOWFL compares these values to the X-Y locations of 
the flow line points in determining the relative position of 
the flow line with respect to any given planform. HEAPSRT 
uses a Heap Sort alogorithm, which is described in reference 

5. 

4. SIFT - Called by HEAPSRT to swap values in the sorting 
process. 

5. CDICLS - Computes the far-field induced drag for simple 
configurations (wing-body) with no dihedral. 
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3.2.4 NEAR-FIELD INDUCED DRAG SOLUTIONS 


This node consists of the single routine, CDRAGNF, (called by WINGAL) 
that computes the near-field induced drag values. CDRAGNF uses 
routine FTLUP in determining the near-field chord force properties. 

3.2.5 SIDE-EDGE SUCTION ANALYSIS 


This node consists of two routines, TIPSUCT and WRTANS. The principal 
routine, TIPSUCT, called by WINGAL, computes the side-edge force and 
values for KV,se for each planform. WRTANS is then called by TIPSUCT 
to compute the values of Kp and KV ,le for each planform. If the 
configuration is cambered/ twisted so that the option for vortex flow 
computation on warped wings is exercised, WRTANS is not called and Kp 
and KV ,le are not calculated. In this situation, TIPSUCT makes 
numerous calls to the LRCGOSF routines for producing the graphics 
output. This section of TIPSUCT will have to be revised considerably 
by users if the VLM program is installed on any other computer. The 
LRCGOSF routines used by VLM are listed in the Appendix A. 
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3.2.6 VORTEX FLOW ANALYSIS 


This node consists of the single routine, VORTEX, called by WINGAL, to 
perform the vortex flow analysis for cambered configurations. VORTEX 
also uses routine FTLUP. 


3.2.7 LONGITUDINAL LOAD DISTRIBUTION 


This node consists of 18 routines, the principal one being CNLONG, that 
called by WINGAL. These routines compute the longitudinal load 
distribution of the configuration. Most of them are associated with 
the Delta Cp or Net Cp surface interpolation from constant Y to 
constant X values. They are (by name only) as follows: 


INTERP 

SUTS 

IQHSCV 

CURV1N 

IQHSD 

CURV12 

IQHSE 

CURVI 

IQHSF 

CEEZ 

IQHSG 

TERMS 

IQHSH 

SNHCSH 

UERTST 

INTRVL 

UGETIO 



The routines beginning with the letters "IQ 11 provide smooth curve 

fitting over randomly distributed data points. IQHSCV calculates the 
interpolating function that is a fifth degree polynomial, and is 
continuous and has continuous first order partial derivatives. This 
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technique is detailed in reference 6. Routines IQHSD through UGETIO 
are support routines for this process. Routines SUTS through INTRVL 
are used to evaluate the the integral of the spline approximation to a 
function of a single variable. These routines were adapted from 
current routines on the NASA — Langley Research Center FORTRAN Math 
Library which is documented in reference 7. This documentation is 
available from the NASA - LaRC Analysis and Computation Division (ACD) 
User Support Office. 

3.3 FILES 


There are nine files associated with the VLM program and these are as 
follows: 

INPUT TAPE6 TAPE20 
OUTPUT TAPE 10 TAPE30 
TAPE 5 TAPE 11 TAPE81 

TAPE81 is equated to INPUT and TAPE 6 to OUTPUT, and all formatted 
INPUT/OUTPUT is performed on these files. Routine READIN reads the 
input from TAPE 81 and prints all of it except the title card over onto 
TAPE5. For the remainder of the program, TAPE5, serves as the primary 
input file. The remaining files are scratch files used by the program 
for auxiliary storage. TAPE 10 and TAPE 11 are used in the matrix solver 
routines and are treated as random access backing store. The CDC 
routines READMS , WRITEMS, OPENMS, and CLOSEMS provide this capability. 
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TAPE10 1 8 reused as a scratch file in the vortex analysis routine. 
There it is employed along with TAPE20 and TAPE30, and subsequently, 
TAPE30 is used to pass information from routine VORTEX back to routine 
TIPSUCT for the final coefficient reporting. 


4 LOAD STRUCTURE 


The tree diagram shown earlier can now be represented as follows: 
GEOMTRY MATXSOL AERODYN CDRAGNF TIPSUCT VORTEX CNLONG 


PROGRAM 

W1NGAL 

where the nodes are main routines Instead of functional descriptors. 
On CDC equipment, this method of loading is called Segmentation and is 
done by the CDC SEGRES Loader Program. The SEGRES program reads a 
directives file, shown in figure 1, and in effect structures the 
relocatable binary file into a set of movable pieces (here, each branch 
of the tree). When the main routine in these branches is called by 
W1NGAL, the Loader brings into memory the entire branch, leaving the 

others out on mass storage. When a branch finishes processing and 
returns control back to the main program (WINGAL), that branch is moved 
back out to mass storage, and another one brought in and loaded into 
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the same space occupied by the previous one. Hence, the same space in 
memory is reused, which reduces the core requirements of the program. 
Complete details regarding the syntax of these directives and the 
actions of the SEGRES Loader Program are given in reference 8. 
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APPENDIX A 

LRCGOSF ROUTINES USED BY VLM 


PSEUDO 

CALPLT 

INFOPLT 

NOTATE 

PNTPLT 


Documentation on these routines is given in reference 9, which is 
available from the NASA — LaRC ACD User Support Office. 
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LISTING OF VLM SOURCE CODE 
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100 FORMAT (1H1»30X*11HINPUT DATA ,H ) READIN21 

200 FOR MAT ( 8 A 10 ) REA0IN22 

300 F0RMAT(I10»1H.»1X>8A10> REA0IN23 

RETURN READ IN24 
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C DETERMINE THE TRAILING EDGE INTERSECTION GE0MT289 

C BETWEEN FIXED AND VARIABLE SWEEP WING SECTIONS GE0MT290 

330 IF ( A8SIRSAR (KJ-SAR (K-l) ) .LT. ( .1/RAO ) ) GO TO 360 GE0MT291 

Y(K)*XPEF(K«-l>-XtK-l)-A(K)*YREF(K+l>+A(K-l)*Y(K-l) GE0MT292 
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SELECT MAXIMUM 8/2 AS THE WING SEMISPAN. IF BOTH FIRST AND GEOMT386 

SECOND PIANFORMS HAVE SAME SEMISPAN THEN THE SECOND PLANFORM IS GEOMT387 
TAKEN TO BE THE WING. GEOMT388 

GEOMT389 

KBOT-1 GE0MT3 90 
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ROOT TREE WINGAL- ( GEOMTRY , MATXSOL , AERODYN , CDRAGNF , TIPSUCT , VORTEX , CNLONG 

,) 


A 

INCLUDE 

GLOBAL 

WINGAL , INFSUB , LOADING , FTLUP , READIN 

ALL , TO THREE , THREFOR,ONETHRE .MAIN ONE , CCRRDD , INSUB23 

* 

GEOMTRY 

a 

INCLUDE 

GEOMTRY, PLANPLT 

MATXSOL 

* 

INCLUDE 

MATXSOL , GIVENS , BLOCKR , TRIANG , SOLVER, BUFFIN 

AERODYN 

a 

INCLUDE 

AERODYN , FLOWFL , CDICLS , HEAPSRT, SIFT 

CDRAGNF 

* 

INCLUDE 

CDRAGNF 

TIPSUCT 

A 

INCLUDE 

TIPSUCT, WRTANS 

VORTEX 

A 

INCLUDE 

VORTEX 

CNLONG 

CNLONG 

CNLONG 

A 

INCLUDE 

INCLUDE 

INCLUDE 

CNLONG , INTERP , IQHSCV , IQHSD , IQHSE , IQHSF , IQHSG 
IQHSH , UERTST , UGETIO , SUTS , CURV IN , CURV1 2 , CURVI 
CEEZ , TERMS , SNHCSH , INTRVL 

A 

END 



FIGURE 1 - SEGMENTATION DIRECTIVES 
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